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GRAVITATIONAL POTENTIAL ENERGY OF THE EARTH: 


A SPHERICAL HARMONIC APPROACH 
Ztevid I%rry Rubincam 

ABSTRACT 

A spherical harmonic equation for the gravitational potential energy of the 
earth is derived for an arbitrary density distriimtion by conceptually brining 
in mass-elements from infinity and building up the earth shell upon spherical 
shell. The zeroth degree term in the spherical harmonic equation agrees with 
the usual expression for the ener^ of a radial density distribution. The second 
degree terms give a maximum nonhydrostatic energy in the mantle and crust of 
-2.77 X 10*® ergs» an order of magnitude below McKenzie's (1966) estimate. 

This figure is almost identical with Kaula's (1963) estimate of the minimum 
shear strain energy in the mantle, a not unexpected result on tile basis of the 
virial theorem. If the earth is assumed to be a homogeneous viscous oblate 
spheroid relaxing to an equilibrium shape, then a lower limit to the mantle 
viscosity of 1.3 x 10*® poises is found by assuming the total geothermal flux is 
due to viscous dissipation. This number is almost six orders of magnitude be- 
low hfocDonald's (1966) estimate of the viscosity and removes his objection to 
convection. If tiie nonequilibrium figure is dynamically maintained by the earth 
acting as a heat engine at one per cent efficiency, then the viscosity is 10** poises, 
a mimber preferred by some (e.g. Cathles (1975)) as the viscosity of the mantle. 
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GRAVITATIONAL POTENTIAL ENERGY OF THE EARTH; 


A SPHERICAL HARMONIC APPROACH 


INTRODUCTION 

Many quantities of geophysical interest, such as the earth's density distri- 
bution, gravitattonai potential, and moments of Inertia, may express^ in 
terms of spherical harmonics or integrals of spherical harmonic coefficients 
(Kaula, 1968). I^nce the earth's gravitational potential ener^ has received 
attention (atul spherical harmonic treatments) in the ^st, particularly from 
MacDonald (1966), McKenzie (1966), and Kaula (1967), it shmild prove worthwhile 
to (terive a ^neral equation for the energy from the viewpoint of spherical 
harmonics. 

We will derive here an equation for the gravitational potential energy for 
an arbitrary density distribution and make a few simple ai^lications with some 
remarks. I^cifically, we will investigate the energy released when a homo- 
geneous earth differentiates into a mantle and core; compute the iK>nhydrostatic 
part of the energy contain^ in the gravity anomalies in the mantle and ci <Mt; 
show its relation to the elastic energy; and estimate a lower limit on the vis- 
cosity of the mantle. We also hope to clear up at least some of the questions 
surrounding the subject of the earth's gravitational potential energy. 
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DERIVATION 


The gravitational potential energy is given 


U = -1g r f dvdv*. 

•'v •'v |r-7*| 

(Sto Blac Donald (1966, p. 230), or an elementary texttxKik on electricity am! 
magiMtism for ^ analogous electrostaUc energy. Tl» p^avitational poteitfial 
energy is seldom discussed in geophysics textbooks.) 

I%e singularity which occxirs when ? « F a{^>ears to make the evaluation 
of Uie integral difficult. We can get around teis pzt^lem, however, in the follow 
ing mani^r. We can conceptually assemble the earUi shell upon spherical shell 
by bringing in matter from infinity and depositing it on Uie earth's surface, 
computing the work necessary to bring in each shell. The sum of the work done 
bringing in all the shells then gives us - U. 

Let us assume for simplicity that the earth is si^erical and nonrotating. 
Neither of these assumptions is restrictive and both will be discuss^ later. 

Let us further assume that the earth's density distrilxition may be expressed 
in terms of normalized spherical harmonics: 


P(S) = /O(r,0, \) = 



(«*) Y^. (0. 


(I) 


with 


(*• = h. (cos0) cos mk, Y^^^ (0, X) = (cos 0) sin mX, 


2 


and (r) « (r), ^ ^^® associated Legei^re 

polyiK>mials, normalized so that 



(cos <^)3 * sin 0 d<^ = 4 - 2 5 q^ 


and 

[Y^ («^. \)]2 dA = 4w 

unit 

sphere 

where dA represents an element of area. 

If we bring in a mass dm = p (r, 0 , X )dAdr from infinity and place it on the 
surface of a partially-assembled earth of mass M and radius r (see Figure 1), 
then the work involved in doing this is V(r, 0 , X ) dm, where V is the gravitational 
|K>tential. Bringing in more masses until we have built a spherical shell of thick- 
ness dr requires work 

-dU = dr J V(r,0. X)p(r, 0, X)dA (2) 

■ur f see 

where dA » r^sin 0 d0dX. (briefly speaking, the layer does work on itself; 
but this is of order (dr)^ in the infinitesimals and may be neglected.) 

The potential V may be expressed in terms of spherical harmonics (Kaula, 
1968, p. 64): 


V(i,«,X) = ™ 


Z] (if 


Yj> . (0.X) 


3 


where A is radial distance to some external poiiU and where 

4tt f x'^2 (x)dx 
Qt -_Zl 

with and The refer of course to the partially- 

assembled earth and sh<xild not be confused wiUi the present-day potential 
confidents. 

The potenUal at A » r may tiien be written 
V (r, 0, \) 2 4t7 G 

'tmi 

aibstituting this expression along with (1) into (2) and working out the integral 
yields 


£o 

(2'C + l) r'^*^ 




-dU= 167 t 2 Gdr 




(2-t + l) 


f 

Ja 


‘'^*’^•{.1 <*>*** 


by the ordiogonaUty of spherical harmonics. 

Adding on more shells until we reach the final radius Rg of the earth gives 
us our expression for the gravitational potential energy; 


'tmi 


where 




( 3 ) 
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We note in passing that our expression differs from McKenzie's (1966) 
equation (10); his expression deals with a layer of mass added to the surface 
of a deformed earth with uniform density. (See Appendix 1). 

Our equation differs also from that of Kaula’s (1967) eq. (3) for the "energy", 
or degree variance, which may be written 

m mi 

The two are related through our eq. (6) (see below). 

RADIAL DENSITY DISTRIBUTION 

As an example of eq. (3), consider the case of a spherical earth with a 
radially symmetric density distribution. Then (r) = (r), = 0 for 

^ 001, and the energy is simply 

U = -16tt^G j PQ(r)r f pg(x)x^dxdr, (4) 

JQ Jo 

where a_ is the mean radius of the earth. (We keep a_ and the final radius R_ 
distinct for reasons made clear below.) This expression agrees with that of 
Urey (1952, p. 174). For an earth with constant density p^ the above equation 
works out to be 

t gm; 

U = - -r = - 2.24 X 10^^ ergs 

where is the mass of the earth and we have used the numerical values given 
in Stacey (1969, Appendix E). 
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An eftrth differentiated into a mantle and cK>re with c^atant respective 
densities and and no change in mass or radius is a little more complex. 


In this case the interior integral in eq. <4) is 


P^ -y for 0 i r i 


and 


aj (r® - a| ) 


Pc 3 +P| 


for a^, < r S 


where a. is the radius of the core. Using this information the energy is 


3 r[n>i“>n2] [2mj-3m2] [2mg + 5 (nij -nij)] 

u = --G.! . 


= - 2.45 X 10^* ergs 


where 


4t7 3 

">1 = — ®c 


4tt 

T 


"*2 = — 


and our numerical values have come once again from Stacey (1969). Hence 
If a hom<^eneous earth differentiates Into the present-day mantle and core, 
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approximately 2. 1 x 10^ ergs is released as heat, in rough agreement with 
Urey (1952, pp. 174-176), who used a somewhat different model of density 
and took compressibility into account. 


ENERGY IN THE MANTLE AND CRUST 

We now turn our attention tc^ the gravitational potential energy; contained in 
the gravity anomalies. 

The density variations po (r) give rise to the observed C of the earth’s 

"vini 'tmi 

gravitational field. By choosing various models for the p a (r) and using the 
known Cj .as constraints, we may use eq. (3) to estimate the gravitational 

■wini 

potential energy in the mantle and crust, where most if not all the gravity 
anomalies are believed to reside. In particular, if we choose p. (r)a (r/R_)”, 
where n is an adjustable parameter, then 


(■•) = <2^ + 1) (n + 't + 3)C.j,^^ r" 


M. 


4t7 r^+3 


as may be found by substituting equation for C^^,; giving 


U 




= _ liiliHllllli) n 

(2n + 5) 


'^mi Rg 


(5) 


If we require that n -• ao , then we are dealing with a surface density distrilxition 
(Dirac delta function), and eq. (5) becomes 


U^* = - C j , (6) 

'val 2 'vmi Rg 
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which is the maximum energy contamed in each harnK>nic» since the material 
is as far from the center of the earth as it can be. The above equation may be 
shown to agree with McKenzie's <1966) eq. (10) after some corrections have 
been made to his expression (see below). 

Most of the nonradial gravitational potential energy resides in the 
term, due to the rrtational equilibrium flattening of the earth. What we would 
really like to know at this point is the energy due to nonhydrostatic equilibrium. 
To find it let uo do the following: we stop the rotation of the earth and assume 
that it relaxes to a spherical shape with radius a^, but retains the same non^ 
hydrostatic C^. as it did before. (This will not be strictly true, of course, but 
makes for a simple case to analyze.) Then the gravitational potential eneigy 
contained in the important second degree {i = 2) harmonics is from eq. (6) 
(going back to the C and S form for the potential coefticients) 

^ ®E 

where we have been careful to use A , the nonhydrostatic part of . Both 
and have been set equal to zero, in accordance with the small amplitude 
of the Chandler wobble. Using the entries in Table 1 of Kaula (1967) as values 
accurate enough for our purposes, namely 

ACjo = 4.70 X 10“* 

= 1.34 X 10'« 

= 2. 40x 10"^ 
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we obtain U*“* = -2.77 x ergs as the maximum amount of ener^ contained 
in the second degree harmonics. This is smaller by an order of magnitude th.an 
McKenzie's (1966) estimate. The discrepancy will be commented upon beiow» 
Probably a reasonable lower limit on the second degree energy can be ob- 
tained from eq. (5) by setting n « 0, so that the anomalous density distribution 
is spread throughout the earth; then 

- 5.44 X 10^® ergs ^ U2^“2.77x 10^® ergs , 

and a guess of ^ -4 x lO*’ ergs is almost certainly right to within a factor 
of 2. 

An estimate of the total gravitational energy in the earth for € 2 can be 

found from Kaula's rule-of-thumb, as given in Kaula (1968, p. 77): 


From (6) 


1 

(2-f + l) 




The series appearing on the right side may be evaluated with the help ot Jolley 
(1961, pp. 64-65, 240), and we have U*** = -6.5 x 10^® ergs. A guess of = 
-1 X 10*® ergs is probably good to within a factor of 2. 
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Oar value of U®** is almost identical with Kaula's (1963) estimate of 
2.94 X 10^^ ergs for the minimum second degree elastic shear strain energy of 
the mantle. This is perhaps not a coincidence; rough equality is expected on 
the basis of the virial theorem, if gravity and elasticity are the two sources of 
potential energy and the velocity of the particles making up the earth aie small 
(see Appendix 2). 

ENERGY OF AN OBLATE SPHEROID 

We have assumed in the derivation for the gravitational potential energy 
that the eartii is spherical and nonrotating. Neither of these conditions hold 
for the real earth, of course. However, the assumptions are not restrictive. 

As far as rotation is concerned, the gravitational potential energy depends only 
on the rela^ve positions of the particles composing the earth and not their 
velocities, hence the rotation of the earth plays no part in the computation of 
the potential energy. Rotation is important, of course, in computing the total 
mechanical energy E of the earth, which is the sum of the kinetic and potential 
energies. For simple rotation about the polar axis this is merely 

E =1 Co;2 + u 

where C is the polar moment of inertia, a; is the angular speed, and other forms 
of energy are ignored. As for sphericity, we may take the earth to be spherical 
by letting R^ be the distance from the center of the earth to the highest point on 
the planet (see Figure 2). In practice this will be the equatorial radius of the 
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earth; thus our distinction between and a^, the two symbols commonly used 
to denote the equatorial and mean radius, respectively. The density of the 
imagined s{^erical earth now happens to be zero in the space between the actual 
surface of the earth and the sphere. This poses no particular problems and is 
no obstacle to a spherical harmonic expansion of the density distribution. 

Let us compute the gravitational potential energy to second degree for a 
homogeneous oblate spheroid as an illustration of this point. The result will be 
used later to estimate a lower limit on the viscosity of the mantle. 

Let the equatorial radius of the spheroid be Rg and the eccentricity be e. 

It may be shown from 





p(r,0, X) Y^. { 4 >. X) dA 


that 


^001 ~ 
?201 (O = ® 


for OS r SRj(l-e2)>^2 


and 


^001 i r 


\ 


for Rg (l-e 2 )<riRg 


^201 “ 2 


(1 


_e 2 ) 3/2 (Rj-r 2 ) 3/2 (R|-r 


2 _ t-2\1/2 


/ 
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where R|^(l - is the polar radius. The other (r) through degree 2 

are zero by the symmetry of the spheroid. 

The above expressions may be substituted in eq. (3) to give and . 
The calculations are tedious in the extreme; only the final results will be given 
here. They are 


U 


boi = 




II 167T^ ^ _2 d5 

*^201 = - GPg Kg 




so that 




( 7 ) 


to the fourth power of e. 

This result is in complete agreement with the exact expression for a homo- 
geneous spheroid given by Lyttleton (1953, p. 36), which is obtained by well- 
known Integral techniques dealing with rotating liquids; 

U = — G p ‘ R; i Arc sin e 

15 E E g 


167t2 ^_2«s r. So 11 4 

= - Gp'R£ 0- — e^- 

15 \ 6 120 
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VISCOSITY OF THE MANTLE 


The shape of the earth is approximately that of an oblate spheroid. The 
flattening factor f, related to the eccentricity through the equation 1 - f = 

(1 - , which best fits the earth in the least squares sense, is f * 1/298.255 

(e.g. Kahn and O'Keefe 1974). If the earth were in hydrostatic equilibrium the 
flattening would be fj^ * 1/299.75 (Kahn and O’Keefe, 1974). Hence the earth is 
flatter than predicted from hydrostatic theory, i.e. the equatorial bulge is too 
big. 

Various mechanisms have been suggested for producing the excess bulge. 
Munk and MacDonald (1960) and MacDonald (1966) thought the excess flattening 
might be a fossil fcailge left over from the remote past when the earth was rotating 
faster. This Implies the earth has a "long memory" (roughly 10^ years), or a 
high viscosity (about 10^® poises), if a linearly viscous fluid is assumed to be 
the appropriate rheology. Goldreich and Toomre (1969) strongly indicated that 
there is no fossil bulge at all; subtraction of the hydrostatic bulge shows the 
earth to be a distinctly triaxial object. They felt that the irregularities in the 
earth’s gravity field might be a by-product of mantle convection and that the 
viscosity was several orders of magnitude smaller than 10^® poises. The 
irregularities would then steer the rotation axis to a position which maximizes 
the polar moment of inertia, thus producing the excess bulge. Wang (1966) 
thought the excess flattening mi^t be due to heavy glaciation at the poles, 
which would squeeze out a bulge due to the weight of the ice. This view was 
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criticized by McKenzie (1966), Kauia (1967), and O'Connell (1971); they indicated 
the ice caps would have to be unacceptably large to produce the extra flattening. 
Kahn and O'Keefe (1974), however, showed that glaciation in Antarctica probid>ly 
produced the gravity field's large third harmonic, giving die earth its "pear 
shape. " Jeffreys (1970, pp. 429-432) felt the earth has finite strength. 

We will not speculate upon the cause of the excess flattening here, but merely 
note that it exists and use it to estimate a lower limit for the viscosity of the 
mantle. Our argument depends upon the excess flattening, the mechanical energy, 
viscosity, and heat flux from the earth. 

Take the earth to be a homogeneous, oblate spheroid with constant density 
and viscosity; thus we will make no distinction between upper and lower mantle 
viscosities. The viscosity may, in fact, be relatively constant throughout the 
mantle (Cathles, 1975, p. 3). 

Let us first assume that the earth has been squashed past its equilibrium 
flattening and is now relaxing back to its equilibrium shape. Heat will be 
generated as the excess bulge subsides through viscous dissipation, subtracting 
energy from E, the total mechanical energy. Hence if we find E, the rate of 
change of mechanical energy, we will have the heat flux due to viscous dissipa- 
tion, by conservation of energy. Let us proceed to do this. 

Let e. , C , ctj. , and denote the equilibrium values of the eccentricy, 
polar moment of inertia, rotational speed, and equatorial radius of the earth, 
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respectively. Let e, C, oj , and refer to the same quantities at some time t. 
Assuming conservation of mass (and volume) 






4t7 _ 


2\l/2 


conservation of angular momentum 


= "h ■ 

and the relationship between the equilibrium eccentricity and angular speed 
(Lyttleton, 1953, p. 38) 


oot 


2nGp. 


3-2e2\ (l-e2) 

(1 arc sin e - 3 


and using eq. (7), we have for the kinetic energy T and gravitational potential 
energy U at time t 


T = 



■ 15 


o;5>b; 





2 p5 
E 


1 



to order e^ , where e Is the only quantity in the above equations which varies 
with time. Adding the two equations together gives 


15 


E = 


167 4 1 /9 9 9^ 

*" + 45 "'"^ J 



^ 16 7T* 

"" “is“ 


G?iRS 


[- ‘ + f^" + !^f 2 +:r <'-*")'] 


(8) 


where we have used 2f = e* and 2f. - e* . 

h ■* h 

Darwin (1879) worked out the relation between the viscosity and the flatten- 
ing. It is for second degree terms 


f=f,+Af 

where 


A f = A f 0 e ^ 


A being the excess flattening at time t = 0, g the gravitational acceleration at 
the earth’s surface, and rj the viscosity. Substitution in eq. (8) yields 


E y 




15 


C Pi «h 




f 

" + 105 



Differentiation of this expression with respect to time yields 


IS 


(9) 


E~ - 


16tt2 


15 




gA 


E 


45 I9ri 


(Af)= 


Taking t to be the present time and Af =* f - f ^ = 1/298. 255 - 1/299. 75 gives 
E y 4.05 X 10^°/t 7 ergs/sec as the present-day heat flow due to viscous dis- 
sipation. The values for f and refer, of course, to the real earth and not 
homogeneous spheroids; but they should be good enough for our purposes. 

Now E must certainly be less than the observed total geothermal flux from 
the earth of 3.15 x 10^° ergs/sec (^acey, 1969, p. 280), which is believed to be 
primarily due to radioactive heating. Thus 


E=4.05x 10 '*Vt 7 e rgs/sec ^ 3. 15 X 10^® ergs/sec 

or 17 > 1.3 X 10^° poises. This is certainly in agreement with the observed 
Fennoscandlan uplif ,, which gives i? = 10^^ poises, the number preferred by 
some as the viscosity of the mantle (e.g. Cathles, 1975, pp. 1-4). 

Instead of relaxing, the excess bulge might be dynamically maintained through 
the earth's action as a heat engine. Goldreich and Toomre (1969) suggested that 
this is indeed the case: the gravity anomalies are a by-product of convection. 
Stacey (1967; 1969, pp. 209-210) estimated the efficiency of the earth's heat 
engine at less than 10 per cent. If one per cent of the earth's heat flux is used 
to maintain the excess bulge, then the viscosity is 10^^ poises, in agreement 
with Cathles (1975). Smaller efficiencies of course yield higher viscosities. 
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It is worth mentioning that Paddack's (1967) observation of the nodal 
acceleration of an earth satellite also permits a lower limit to be put on the 
viscosity. The nodal rate of i satellite is proportional to = -Cjq term 
in the gravity field expansion. Since is also proportional to f for a homo- 
geneous oblate spheroid, we have = - Kf, where K is a constant of no concern 
to us here. Differentiation with respect to time and division by Q, yields 

_Af 

n ■ f i9t) • 

s 

Paddack found the left side to be zero to one part in a million. This gives 
7? > 2 X 10*® poises, which is smaller than our limit by five orders of magnitude. 
Thus our heat flow argument puts a more stringent lower limit on the viscosity 
than does the observation of the satellite. 

DISCUSSION 

McKenzie (1966) estimated the second degree gravitational potential energy 
at -2.02 X 10^° ergs, a factor of ten larger than our estimate of -2.77 x 10^^ 
ergs. We can clear up this discrepancy by noting three things about McKenzie's 
analysis. First, there is a spurious factor of 4 t 7 in his equation for dm (the 
third equation above his eq. (10)). Second, his gravitational potential coefficients 
C-Ci ^**®^**S °***^ notation) differ from our normalized coefficients by the 
relation = (2 - )*'^^ C^. . He erred in using the numerical vahies of 

for his values of C^. in computing the energy. Last, his procedure for 
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deriving the gravitational energy is slightly inconsistent: in some places in the 
derivation a spherical surface is assumed and in others a deformed surface. 

If a spherical surface is consistently assumed (i.e. using his eq. (6) for the 
potential coefficient), and the above factors of and (2 - m-e noted, 

then his eq. (10) becomes identical with our eq. (6). 

MacDonald (1966) after a lengthy analysis estimated that 2 x 10^"^ ergs of 
gravitational potential energy was pent up in the earth, a factor of 10® larger 
than our estimate of total energy release of 2 x lo^® ergs, which may be derived 
from our eq. (9). Use of MacDonald’s figure in our viscosity argument results 
in 7) = 10^® poises, in agreement with the number given on his page 227. This 
assumes, however, that the entire geothermal flux comes from viscous dis- 
sipation. Attributing part of the heat flow to radioactive heating results in even 
higher values. Such viscosities rule out convection on a meaningful Ume-scale; 
the mantle would not overturn even once in the entire history of the earth. 
Adoption of our much lower number does not necessarily rule out convection. 

But we can make no statement as to the actual value of the viscosity: in principle 
it could be any value above 1.3 x 10^° poises. However, a thermodynainic 
efficiency of one per cent yielding 10^ ^ poises is suggestive. 

We can answer two questions raised by Kaula (1967, pp. 790 and 792). 

First, he wondered whether the energy associated with the gravity anomalies 
would become available as heat should the anomalies disappear. The answer 
depends delicately on the shape of the earth and its internal density distribution. 
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The energy associated with the relaxation of an oblate homogeneous spheroid 
produces heat, as we have found in the last section; i.e. the mechanical energy 
decreases. But the gravitational potential energy increases when the anomalies 
disai^ear in a spherical earth, as evidenced 1^ our eq. (5). Therefore computa- 
tions involving the gravitational energy must be done with great care. Second, 
Kaula wondered whether the harmonics in a shell of matter interacted only with 
themselves to produce self-energy. Oir derivation answers this question in the 
negative: the gravitational energy arises through the interaction of each har- 
monic of a shell with the same harmonic in the other material composing the 
earth. 

We hope to hai'e answered at least some of the questions surrounding the 
topic of the earth's gravitational potential energy. Future investigations should 
prove fruitful in ciarifying still further our understanding of the physics of the 
earth. 
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Figure 1. Bringing in a mass ^ from infinity to the 
surface of the partially-assembled earth of mass 
M and radius r. 
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FIGURE CAPTIONS 


Figure 1. Bringii^ in a mass from infinity to the surface of the partially 
assembled earth of mass M and radius r. 

Figure 2. Imaginary sphere containing tbe earth. Rg is the distance from the 
center to the hipest point on the earth. 


APPENDIX 1 


McKenzie's {1966) equation (10) is, in his own notation, 



( 28 + 1 ) Cjo 

t 


+ 22 (*-l) (28+1) 

8,n»#o 



where g is the gravitational acceleration at the surface of the earth, and a 
and M are the mean radius and mass of the earth, respectively. The 
and are his potential coefficients. 

The equation can be made to ^ree with our eq. (6), after corrections, as 
stated in the text. Our eq. (6) can also be derived from the equation for the 
potential energy for a surface distribution: 


.-if 


aV dA. 


Here a Is the surface mass density and V is the potential. 


APPENDIX 2 


The virial theorem states (Goldstein, 1950, p. 70): 

-2 <T) = 

i ^ 

where T is the kinetic energy of the earth, F| and are the total force on 
and the position of the ith particle of the earth, respectively, and the angular 
brackets denote time - aver^es. We will henceforth drop the brackets. 

For small perturbations in the positions of the particles the above equation 
becomes 

-2 5T =]^ Fj • 
i 

= 5^?!^ • 6?j+5]Fj^ • 6?, 

i i 

where we have decomposed the total force on each particle into the gravitational 
and elastic forces. The right side of the equation is nothing more than the work 
done, which is equal to minus the change in energies, so that 


2 fiT = 6U + 6W, 


where U is the gravitational energy and JW is the elastic energy. 
If 


6T 


« su 


A2-1 


then 


6U as - 5 W, 


which is the statement desired. 

To test the inequality, we note that 


6U 




^ 3 X 10^ ergs 


For 6T to approach jU 2 ”“’‘l in value, we would have to have particle speeds 
V such that 


6T 25 



2s 3 X 10® 


ergs 


or V = 10 cm/sec. Since V is more like 2 cm/year, we can assert the in- 
equality with confidence (assuming heat plays no role in the energy balance). 
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